Transverse mass dependence of two-pion correlations 
in Au+Au collisions at ^snn = 130 GeV 
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Two-pion correlations in -^snn = 130 GeV Au+Au collisions at RHIC have been measured 
over a broad range of pair transverse momentum kr by the PHENIX experiment at RHIC. The 
kr dependent transverse radii are similar to results from heavy ion collisions at ^snn = 4.1, 4.9, 
and 17.3 GeV, whereas the longitudinal radius increases monotonically with beam energy. The ratio 
of the outwards to sidewards transverse radii (-Rout/-Rsidc) is consistent with unity and independent 
of kr- 

PACS numbers: 25.75.Dw 



The influence of Bose-Einstein statistics on the corre- 
lation of identical charged pions at low relative momen- 
tum was first used to probe the space-time structure of 
pion emission in pp annihilations and has subsequently 
been applied to relativistic heavy-ion collisions from the 
Bevalac to RHIC (see for a recent review), and to 
a wide range of systems including e^e~ annihilations 
The correlation function is defined as the ratio of the 
two-particle probability distribution to the product of the 
single particle distributions. For a static source with no 
final state interactions, it is related to the Fourier trans- 
form with respect to q = pi — P2 of the source distribu- 
tionp(r), P(pi,p2)/P(pi)P(p2) = l + |p(q)P §■ If the 



source is parameterized as a multi-dimensional Gaussian, 
the enhancement in the correlation function is a Gaus- 
sian, and the Gaussian widths are each inversely propor- 
tional to the source dimensions in the canonically conju- 
gate spatial variables. The extracted source dimensions 
are commonly referred to as HBT radii, after a simi- 
lar technique developed by Hanbury-Brown and Twiss 
to measure stellar radii [|o|. For dynamic sources, such 
as rapidly expanding sources in heavy-ion collisions, the 
correlation function measures "lengths of homogeneity" , 
or the relative separations of the pions with low rela- 
tive momentum. This leads to source radii which depend 
strongly on kx, the mean transverse momentum of the 
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pion pair ||Tl|-|l^. If the dynamics are correctly mod- 
eled, then both the source geometry and rate of expan- 
sion can be deduced by measuring the kr dependence of 
the radii. The existence of a connection between HBT 
radii and heavy-ion source geometry is established by 
the dependence of the radii on system size cen- 
trality [§,D, and reaction plane Interest in Bose- 
Einstein correlations in heavy-ion collisions is driven by 
the expectation that HBT radii are sensitive to the large 
and/or long-lived sources which may accompany a QCD 
phase transition ■ Recent calculations predict that 

the greatest sensitivity to a long-lived source will come 
from measurements of the correlation function at high 
kr (> 0.3 GeV/c) 

We present new measurements from the PHENIX ex- 
periment on two-pion correlations in Au-|-Au collisions at 
v^siiT = 130 GeV in the region \ri\ < 0.35, 0.2 < kr < 1.0 
GeV/c, significantly extending previous measurements 
by STAR up to a mean-fcr 0.63 GeV/c. The data 
are compared to theoretical predictions for RHIC and to 
HBT radii from lower energy collisions at the SPS and 
AGS. The kx dependence of the transverse radii is used 
to extract a geometric transverse radius. 

The PHENIX experiment has been described in detail 
elsewhere |^,^ . For this analysis we utilize a subset of 
the detectors in PHENIX. We use the hadronic particle 
identification capabilities present in the west arm of the 
PHENIX spectrometer perpendicular to the beam direc- 
tion with polar and azimuthal ranges of \r]\ < .35 
and 7r/4, respectively, during its first year of running. 
In this analysis, the vertex is determined with a zero 
degree calorimeter (ZDC) and a pair of Cerenkov beam- 
beam counters (BBC). Pattern recognition and momen- 
tum reconstruction rely on a drift chamber and a pad 
chamber which occupy the region between 2.0 and 2.5 
meters from the beam axis. The momentum resolution 
from these detectors is Sp/p = 0.6%©3.6%p. Particle ve- 
locity is determined from the differential time measure- 
ments of the BBC and the electromagnetic calorimeter 
(EMC) with a combined rms resolution of 600 ps, 
coupled with the path length determined from pattern 
recognition. The momentum determination and particle 
identification method are similar to |Q , except that the 
time of flight is measured by the EMC. A pion is defined 
as being within 1.5 standard deviations of the pion mass- 
squared peak but at least 2.5 standard deviations away 
from the kaon peak. After applying inter-detector as- 
sociation cuts the background from mis- associated EMC 
hits is ^ 10% as determined by a hit randomization tech- 
nique. This background does not significantly distort 
the extracted radius in the correlation measurements, al- 
though it reduces the measured correlation strength (A). 
We did not correct for this background in our correlation 
analysis. 

A total of 493K events in the most central 30% of the 
cross section survive all offline cuts. This sample con- 



tains 3.1 million tt"*" pairs and 3.3 million tt pairs in the 
analysis, and has a mean centrality of 10%. 

The pion correlation function is determined from pairs 
of identical pions. The normalized probability of detect- 
ing two particles with relative momentum q = Pi — P2 
and average momentum k = (pi -I- p2)/2 is determined 
experimentally by the ratio of pairs from the same event 
(A) with those from different events (B): C2(q, k) = 
y4.(q, k)/i?(q, k). Pairs of particles within 2 cm of each 
other in the drift chamber are eliminated from the anal- 
ysis in both the real and background samples. Pairs that 
share the same EMC cluster are also removed from both 
samples. Finally, all pairs in the mixed background sam- 
ple are required to be from events with a reconstructed 
BBC collision vertex within 1cm of each other. 

We correct for the Coulomb interaction of the pairs in 
the correlation function by parameterizing the source as 
a Gaussian distribution in the pair center-of-mass frame 
and performing an iterative procedure |2^ ] which ac- 
counts for the finite resolution of the detector. This pro- 
cedure applied to the distribution of 7r+-7r^ pairs is in 
agreement with the data, although the statistics in the 
Run-1 opposite-signed analysis are not sufficient to in- 
dependently determine the required Coulomb correction. 
Systematic studies of the Coulomb correction which vary 
both radius and magnitude within reasonable constraints 
produce variations in the final radii which never exceed 
0.25 fm. 

The relative momenta are projected into the variables 
Qiong, along the beam direction, gout, parallel to the trans- 
verse momentum of the pair kr = |(pti + Pt2)i and 
9sidc, perpendicular to giong and Qout |n|,|l|l. These vari- 
ables are calculated in the longitudinal co-moving system 
(LCMS), obtained by a longitudinal boost from the lab 
frame to the frame in which the longitudinal pair veloc- 
ity vanishes. This frame is commonly used for sources 
expected to be invariant under longitudinal boosts [ p6[ . 

The fully corrected correlation function for tt~ pairs 
is shown in the top panels of Fig. Q the large q region 
of the correlation function has been normalized to 1 in 
the plots. The data are fit to a Gaussian parameteriza- 
tion of the source using a MINUIT based log-likelihood 
method 

C2 — 1 + Acxp( — i?iong^'7long^ — -Rsido^gside^ — ^out^gout^) 

(0.1) 

where i?iong) -Rside, and i?out are the conjugate variables 
to giong, gside, and gout, respectively. Errors quoted in the 
tables and figures are statistical only. Systematic errors 
come mainly from the Coulomb correction and depen- 
dence of the results on the two-track distance cuts. The 
combined systematic error for these effects, estimated 
by varying the cuts and corrections within reasonable 
bounds, is 8% for i?iong, and i?sidc, and 4% for i?out- The 
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systematic error from residual correlations in the event- 
mixed background |^ is 2%, yielding a total systematic 
error of ~8% for i?iong and i?side and ~4.5% for i?out- 

The data set is subdivided into three kx bins of equiv- 
alent statistics in order to study the momentum depen- 
dence of the correlation function. In Fig. ^, the radii 
for 7r~ pairs are shown to agree within statistical and 
systematic errors with previous measurements for over- 
lapping kx bins at this energy for the 12% most cen- 
tral events. For STAR, the mean pair centrality can 
be approximated by the geometric mean of 8%, which 
is slightly more central than the mean pair centrality 
of 10% for the PHENIX data. This figure also shows 
kx dependent radii for mid-rapidity pions from central 
collisions for = 17.3 GeV Pb+Pb §,0 and for 

y/s^ = 4.9 and 4.1 GeV Au-hAu ^j^. For the trans- 
verse radii, i?out and i?side, the variation with collision 
energy is generally smaller than the statistical and sys- 
tematic errors of the individual data points. There is 
no evidence for a change in the low-fc^ extrapolation of 
i?sidc with increasing ^snn which would indicate a larger 
geometric source at higher energy. Nor is any change ev- 
ident in i?out relative to -Rsidc at high-fcr, indicating a 
longer-lived source. This result is surprising given the 
factor of change in the total charged particle mul- 
tiplicity per unit rapidity at mid-rapidity Only 
i?iong exhibits a significant variation with collision energy. 
To quantify this difference, we fit the i?iong dependence 
to A/ y/rriT p^ , p^ , p9t for the three sets of beam energies. 
The results arc overlayed with the data in the bottom 
panel of Fig. | and yield A = 3.32 ± 0.03, 3.05 ± 0.06, 
and 2.19 ± 0.05 fm-GeV^ for = 130, 17.3, and 

4.9/4.1 GeV respectively. 

Although a finite emission duration contributes to 
i?out but not to i?sido, dynamical correlations affect the 
two radii differently. A quantitative determination of the 
source lifetime can only be performed in the context of 
a dynamical model. The lower panel of fig. ^ shows the 
kx dependence of the ratio i?out/^sidc for PHENIX and 
STAR along with recent calculations for a thermalized 
source which undergoes a first order phase transition 
at critical temperatures (TJ of 160 and 200 MeV [|o). 
The rise in -Rout/^sidc which comes predominantly from 
a hadronic re-scattering phase is not present in the data, 
and the values of 1.6 (T^ = 160 MeV) and 2.2 (T^ = 
200 MeV) at high kr are excluded. 

An additional consequence of strong dynamics occurs 
for sources in which the transverse expansion is relativis- 
tic. In this case, i?out measured in the LCMS frame is 
Lorentz contracted by the 7 of the pion source velocity 
along the direction of gout [^70- Current Lorentz in- 



variant formulations of the correlation function 14 
are insufficient to determine the source velocity due to 
transverse expansion, however, the pair center-of-mass 
system (PCMS) can be used to provide an upper limit 
on i?out p3]- The correlation function for tt^ pairs in the 



PCMS frame is shown in the bottom panels of Fig. ||, 
and fit results for i?^^^^ are listed in Table |. As ex- 
pected, i?sidc and -Riong are equal to the corresponding 
LCMS parameters within errors. 

Two analytic expressions have been used to describe 

-Rsidc as a function of = 
versely expanding source, 



rni for a trans- 



-Rsidc^(TOT) 



-Rsidc^(™T) 



Rcr 



1 + 



(0.2) 
(0.3) 



Eq. 3.2 is a first order approximation in for a lon- 



gitudinally boost invariant source with finite tempera- 
ture, T, and expansion velocity, /St — /?/ p/flgcom, where 
Rgcom is the Gaussian transverse radius |1J]. Eq. 0^ 
includes an additional term in the approximation and 
the linear transverse expansion velocity is replaced by a 
transverse rapidity, rjT = ?//p/i?gcom fl^- For a trans- 
verse surface rapidity of ry/ = 0.85 (/?/ — 0.69) and 
T = 125 MeV a fit of Eq. lO to the PHENIX 



-Rsidc niT dependence yields, Rg 



.1 ± 0.3 fm with a 



X /dof ~ 9.6/6. To assess systematic errors the PHENIX 
2|, yielding i?gcom = 6.7±0.2 fm 



I • 6.87 = 

5 



data are also fit to Eq 
and xVdof = 9.1/6, andthe STAR data are fit to Eq. 
yielding i?geom = 9.4±0.1 fm with xV^of = 21/6. These 
fits are shown in the top panel of Fig. ^. All values of 
J^gcom are significantly larger than the comparable ID 

rms radius for a Au nucleus p6[ of 
3.07 fm. 

In conclusion, we have extended the measurement of 
two particle correlations for Au+Au collisions at y/sm^— 
130 GeV to < kr > ^ 0.63 GeV/c using the PHENIX 
detector at RHIC. Values of i?^^^^ are used to con- 
strain the Lorentz effects for a relativistic transverse ex- 
pansion. Fitting i?sidc(fcT) to two analytic expressions 
for an expanding source yields a transverse geometric 
radius that is much larger than the comparable radius 
for Au. We find that Riongikr) increases monotonically 
with collision energy, yet no energy dependence is dis- 
cernible in the kx dependence of i?out and -Rsidc, and 
the ratio, -Rout /-Rsidc, is consistent with unity and in- 
dependent of kr- The results for the transverse radii 
are contrary to common expectations for a first order 
phase transition in Au+Au collisions at these energies, 
as demonstrated by the comparison to a typical hydro- 
dynamic model with hadronic-rescattering. Therefore, 
we conclude that current concepts regarding the space- 
time evolution of the pion source inferred from two-pion 
correlations in Au-f Au collisions at RHIC will need to be 
revised. 
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FIG. 1. The three dimensional correlation function for tt^ 
pairs versus giongjQsido, and gout in both the LCMS frame (top) 
and pair center-of-mass frame (bottom). The data are plotted 
versus one momentum difference variable while requiring the 
other two to be less than 40 MeV/c. The lines correspond to 
the fit to the entire distribution. 
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FIG. 2. HBT radii for pion pairs as a function of kr mea- 
sured at mid-rapidity for various energies from E895 (^snn ~ 
4.1 GeV), E866 {^^=4.9 GeV), NA44, WA98 {yd^=17.3 
GeV), STAR, and PHENIX (yi5^^=130 GeV). The bottom 
plot includes fits to Aj ^Jmr for each energy region. The data 
are for tt^ results except for the NA44 results, which are for 
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FIG. 3. The top panel shows the measured i?sidc from iden- 
tical pions for STAR and PHENIX. The line solid line is a fit 
of Eq. ^ to th e P HENIX data, and the dashed line is t he 
same fit for Eq. 0.2. The dot-dashed line is a fit of Eq. D.3 to 
the STAR data. The bottom panel shows the ratio -Rout IRa&c 
as a function of fcr overlayed with theoretical predictions for 
a phase transition for two critical temperatures. 



TABLE I. The dependencies of the tt"*' and tt" radii in 
the LCMS and PCMS frames. All momenta are in MeV and 
all radii are in fm. The errors are statistical only. 





kT (MeV) 


200 - 400 


400 - 550 


550 - 1000 






333 


472 


633 




^inv 


6.74 ±0.31 


6.42 ± 0.46 


3.46 ± 0.46 






0.423 ± 0.037 


0.389 ±0.039 


0.287 ± 0.048 




-Rlong 


6.01 ±0.45 


4.76 ±0.35 


2.97 ±0.38 






4.81 ±0.30 


3.74 ±0.36 


2.79 ±0.37 




Rowt 


4.78 ±0.30 


3.76 ±0.26 


2.59 ±0.46 




pPCMS 
-fl-out 


11.35 ±0.69 


12.20 ± 1.02 


8.60 ± 1.13 






6.00 ± 0.30 


5.96 ±0.41 


4.58 ±0.48 






0.431 ±0.079 


0.405 ±0.067 


0.353 ± 0.062 




-Rlong 


5.69 ±0.76 


4.77 ±0.49 


3.76 ± 0.41 




^?side 


4.67 ±0.38 


4.13 ±0.45 


3.22 ±0.35 




^out 


4.69 ±0.58 


3.75 ± 0.40 


2.81 ± 0.34 




pPCMS 
-fl-out 


11.27 ± 0.72 


12.42 ± 1.18 


11.89 ± 1.73 
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